Extracting spectral density function of a binary composite without a-priori assumption 
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The spectral representation separates the contributions of geometrical arrangement (topology) and in- 
trinsic constituent properties in a composite. The aim of paper is to present a numerical algorithm based on 
the Monte Carlo integration and contrainted-least-squares methods to resolve the spectral density function 
for a given system. The numerical method is verified by comparing the results with those of Maxwell- 
Garnett effective permittivity expression. Later, it is applied to a well-studied rock-and-brine system to 
instruct its utility. The presented method yields significant microstructural information in improving our 
understanding how microstructure influences the macroscopic behaviour of composites without any intri- 
cate mathematics. 
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Theory of mixtures and their electrical properties have 
attracted researchers to seek a relation between intrinsic 
properties of the parts forming the mixture (constituents) 
and their spatial arrangement inside the mixture 
Bergman|4] has proposed a mathematical way for repre- 
senting the effective dielectric permittivity of a binary 
mixture as a function of permittivities of its constituents, 
£i and 82, and an integral equation, which includes the 
geometrical contributions. It is called the spectral den- 
sity representation (SDR). After the introduction of non- 
destructive measurement techniques and systems, such as 
electricalQIsIl or acoustic impedance spectroscopy |3|, the 
impedance of materials (either pure or composite) could be 
recorded for various frequencies u. Then, the frequency 
could be used as a probe to obtain microstructural informa- 
tion with the application of the SDRlHHHHiEO- 
This can only be achieved if (/) no influence of v on the 
geometrical arrangement of phases is present lfl2il . and (//) 
the intrinsic properties of phases are known as a function of 
V. Numeric all 8. A JCM and analytical|,5|, (6|, |7| approaches 
have been used and proposed to resolve the spectral density 
function (SDF) for composites. Although numerical ap- 
proaches could be prefered over the analytical ones, which 
are emprirical expressions and are not universal, they solve 
a nontrivial — ill-posed — inverse problem|9|. Here, we ap- 
ply a recently developed numerical method|13] to extract 
the SDF of a binary mixture. The method is based on the 
Monte Carlo integration and constrained-least-squares (C- 
LSQ) algorithms. By using this procedure the integration 
constant becomes continues rather than discrete as in reg- 
ulation algortihms. First, the verification of the proposed 
method is presented by considering the Maxwell-Garnett 
(MG) effective dielectric functioni 14]. Later, it is applied 
to the dielectric data of a rock-and-brine system| 15], which 
has been also used by Refs^ and to test their analytical 
expressions. 

For a binary composite system with constituent permit- 
tivities £1 and E2, and concentrations qi and q2, (qi + ^2 = 
1), and with an effective permittivity e^, the SDR is ex- 



0.4 




0.2 0.4 0.6 0.8 



FIG. 1 : Parametric plot of the scaled mixture permittivity. The 
symbols are the analytical model of Maxwell-Garnett equa- 
tion, and the solid lines ( ) are the values calculated from 

the spectral functions obtained from the proposed numerical 
method. The semi-circles from large to small corresponds to 
92 ~ {0.95, 0.7, 0.5, 0.3, 0.05}, respectively. The inset is the en- 
largement of the values close to the origin for 92 = {0.05, 0.30}. 

pressed as lfl^ . 

Ae,/A,, -A^= [ g,(x) [1 + er'Aj.x]-' dx (1) 
Jo 

where, = — Sj, and is complex and frequency de- 
pendent. Aj is a constant, and depends on the concentra- 
tion and structure of the composite. The SDF is g{x), and 
it is sought by thepresented procedure. The SDF satisfies 
/ gj{x)dx = qjf^,^ld\ and / xgj{x)dx = qj Qi/d, where 
d is the dimension of the system. The shape of the inclu- 
sions in a matrix can also be related to d lfivl] . Finally, x is 
called the depolarization factor. 

The numerical procedure is briefly as follows: first the 
integral in Eq. (0 is written in a summation form over 
some number of randomly selected (known) x„-values, 
x„ € [0, 1]. This converts the non-linear problem in hand 
to a linear one with gj^ being unknowns. Later, a C-LSQ 
is applied to get the corresponding gj^ -values: 

mill ||A - Kgj„||2 and gj„ > (2) 
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FIG. 2: Calculated spectral density distributions, which cor- 
respond to delta sequences. The spectral functions from left 
to right corresponds to 92 = {0.95, 0.7, 0.5, 0.3, 0.05}, re- 
spectively. The correcponding (calculated) A2 values are 
{0.002, 0.012, 0.029, 0.064, 0.358}, respectively, for the consid- 
ered concentrations. The dashed lines ( ) show the positions 

of the actual delta-functions for the MG expression. 

where A is the left-hand-side of Eq. Q, and K is the 
kernel-matrix, [1 -|-e^^AjiX„]^^. When this minimization 
is run over-and-overwith new sets of x„-values, most prob- 
able gj^ -values are obtained. For a large number of mini- 
mization loop, actually the x-axis becomes continues — the 
Monte Carlo integration hypothesis. Filially, the weighted 
distribution of g^^^ versus x„ leads gfx^ lflal . 

Application of the numerical procedure to the MG 
expression should yield delta function distributions for 
gfxllsL [Toll . The dielectric function for a d-dimensional 
(or composite with arbitrary shaped inclusions) MG com- 
posite is defined as 

Ee = ei[l + dq2 A21 (gi A21 + dsi)-^]. (3) 

The resulting SDF is then, 

g,{x)=6[x-{l-q,)/d]. (4) 

We choose the following values for dielectric functions of 
the phases: Si = I — ^(100eow)"^ and £2 = 10 — 
i{eoUj)~^ with LO = 27ru and Sq = 8.854 pF/m. The 
left-hand-side of Eq. (0 without the constant A2 is plot- 
ted for a 3-dimensional composite (d = 3 which cor- 
responds to spherical inclusions) in Fig. Q as a parame- 
teric plot of the imaginary part of Aei/A2i against its real 
part. The graph is a semi-circle for the MG expression. 
In the figure, five different concentration levels are plotted, 
52 = {0.05, 0.3, 0.5, 0.7, 0.95}, the inset shows the en- 
largement close to the origin, which illustrates the low con- 
centrations, q2 = {-05, 0.3}. The size of the semi-circles 
are proportional to the concentration of Phase 2. The analy- 
ses performed on the scaled effective permittivity, Eq. Q, 
with the help of the applied method yield the solid lines 

( ) in the figure. 

The corresponding g(x) are plotted in Fig. |2lon a log- 
log scale. In the figure, the expected locations of g{x) 

from Eq. Q are also shown with dashed lines ( ). 

The g( a; ) -distributions obtained are analized by the Levy 
distribution! 19], which generates a delta-squence L20J . The 



TABLE I: Compaiison between the reults of the proposed nu- 
merical approach and those of the Levy statistics and the given 
analytical SDF for the MG effective permittivity expressions for 
various concentrations. The bars on the quantities indicate that 
they are calculated from the numerical results. 
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"Calculated using the resulting g2(a;). Known from the definition of 
gj (x) — integral gj (a;)da; is equal to this value. 

''The localization parameter for the calculated Levy distribution. The shape 
parameters and the amplitude of the Levy distributions are disregarded. 

''Known from the definition of the SDF for the MG expression, Eq. PI . 

'j42-value calculated before the numerical procedure using Eq. 0. 

''Mean A2 -value calculated during each Monte Carlo integration step in the 
numerical procedure, Eq. J2|- 

''Calculated using the resulting g2 {x) and x-values. Known from the defi- 
nition of gj (x) — the values is equal to the integral Sxgj {x)dx. 
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FIG. 3: Parametric plot of the scaled rock-and-brine permittiv- 
ity. The symbols (•) are the experimental data of Refil5i The 
chain line (- ■ -) is the results for the same assumptions as Refs0 

and 6. The solid ( ) and dashed ( ) lines are results 

obtained by two different water and composite conductivities; 

(72 = 0.85 S/m and oe = 0.041 S/m ( ), and 02 = 0.85 S/m 

and (Je = 0.038 S/m ( ). 

solid lines ( ) illustrate the appropriate Levy distribu- 
tions. Various parameters from the statistical analyses and 
their expected values are presented in Table |2 The con- 
centration values, q^, calculated from the integration of 
g(x) without a-priori assumption are < 1% for the con- 
sidered higher concentrations, and it is around 5% for the 
lowest concentration, q2 = 0.05. The localization pa- 
rameter for the depolarization factor x, which is the most 
propable depolarization value, can be calculated by the in- 
tegration of [1 — g{x)\/d or with the help of statistical 
analysis. The estimated depolarization factors x are within 
< 1% of the actual values stated by the proposed analyti- 
cal expression|l5lllfl!|. Finally, the product of the concentra- 
tions qiqi, the integration of 2>xg{x), calculated have also 
very good aggreement with those values expected from the 
definitions of the SDR. 
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We also test our procedure on a rock-and-brine com- 
posite svstem Ll51 . which has been studied by various 
scientists pJl 01 . The same assumptions as in Ref|5|and0 
are made to calculate the dielectric function of the brine 
(water-salt solution). The ohmic conductivity of the water 
is taken to be a2 = 0.93 S/m, later the dielectric function 
of the brine at T = 75 °C is calculated by the following 
expression! 5, 6], 



4 m 

S{T) 
e',{T,S) 
e2{T,S) 



94.88 - 0.2317r + 0.000217r^ 
5.363 [(T + 7)(82a2)"^ - 0.0123]-^ "''^ 
[e^Q-^ + 0.04175(1000 - S)-^Y^ 
e',{T,S)-ia2{eoUj)-' (5) 



The relative permittivity of the rock is taken to be con- 
stant without any imaginary part, Ei = 7.5. The result- 
ing scaled dielectric quantity in Eq. Q is presented in 
Fig. |3l Similar to Fig. Q a semi-circle-like shape is ob- 
served. The first analysis with the above considerations 
results in an unsatisfactory calculated as presented with 
the chain line (- • -) in Fig. |3 The low frequency side 
(uj < 30 MHz) of the real permittivity has discrepancies. 
Therefore, the experience of the author regarding dielectric 
data analyses suggests that the measured values at the low 
frequencies do not particularly satisfy the Kramers -Kronig 
relations lfl3l I21I1 . Consequently, the application of the 
Kramers-Kronig relations yield lower effective compos- 
ite conductivity then the original data, = 0.055 S/m. 
Therefore, two different conductivities are adopted = 
0.041 and 0.038 S/m, while we keep the conductivity of 
the water constant and lower than the previous consider- 
ation = 0.85 S/m. With these parameters as inputs, 
the resulting effective permittivity values have better agree- 
ment with those of measurements. And if compared to the 
results of Stroud et al. |5] and Ghosh and Fuchs |6], our 
values have less residual than theirs. 

In Fig.|5l the obtained g{x) are presented. It is striking 
that two very distinct peaks are observed whatever the ini- 
tial assumptions for the conductivities of the water as well 
as the composite are. The g(x) can be divided into three 
sub-SDF, which are located around X = {0, 0.004, 0.04}. 
It is clear that the original data can be modeled by only 
two SDF as delta sequencesi 10, 11, 22] without a sophisti- 
cated mathematics. The SDF of Ref 6 is also displayed as 
a comparison with the thick chain line (- • -), which has 
been valuable to give limits for the depolarization factor 
X. However, in the case of Kenyon's data lil51 it overes- 
timates the upper limit, which has been 1. The two pe- 
culiar depolarization factors resolved from the peaks have 
concentrations of 0.111 and 0.023, respectively, which are 
calculated from the Levy distributions. The low x-side 
of g(x) yields a very small concentration (~ 10~^) for 
that particular depolarization process. If we take into ac- 
count the yielding concentrations of the brine in the system 
(^2 ^ 0.134), we can state that the three peaks correspond 
to oblate to needle like porous structures of the brine with 
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FIG. 4: Measured (• from Ref. 15) and re-calculated dielec- 
tric permittivity K(ee) and alternating current conductivity ae = 
Q'(ee£o'^)- The chain line (- ■ -) is the results for the same as- 
sumptions as Rets 1^ and 1^ The line legends are the same as in 

Fig.m 




FIG. 5: Calculated spectral density distributions. The lines rep- 
resent the appropriate fitted Levy distributions, and their legends 
are the same as in Fig.|3l The thick chain line is the SDF g{x) of 
RefS 



shape factor estimates d ~ Qi/x w |10^, 200, 1.2} for 
X = {0, 0.004, 0.04}, respectivelytl?]. It is clear that the 
brine-phase forms channel-like structures, because of high 
d values. Continuous percolating paths are formed when 
g(x) = 5{x) and d — > 00, which corresponds to a struc- 
ture with channels parallel to field direction. 

As a concluding remark, an effective numerical method 
is presented to extract the SDF of a binary composite sys- 
tem. It is tested on both 'ideal' and measured dielectric 
data for composites. The proposed method not only ex- 
tracts the SDF, it also yields volume fractions of both con- 
stituents as well as the correponding depolarization pro- 
cesses and phase shapes even if they are not known at the 
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beginning. It is shown that it can resolve unique individual 
depolarization processes, which could indeed be used to 
obtain valuable microstructural information regarding the 
composite and its constituents in various research fields, in 
which impedance spectroscopy is used for characterization 
of materials, such as, polymeric, pharmaceutical, biologi- 
cal, building, colloidal, porous, etc. 
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